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ff^ . We study the electron-positron system in a strong magnetic field using the differential Bethe- 

o , 

, Salpeter equation in the ladder approximation. We derive the fully relativistic two-dimensional 

O 

form that the four-dimensional Bethe-Salpeter equation takes in the limit of asymptotically strong 
constant and homogeneous magnetic field. An ultimate value for the magnetic field is determined, 
r~| , which provides the full compensation of the positronium rest mass by the binding energy in the 

maximum symmetry state and vanishing of the energy gap separating the electron-positron sys- 

><' 

■ tem from the vacuum. The compensation becomes possible owing to the falling to the center 

phenomenon that occurs in a strong magnetic field because of the dimensional reduction. The 
solution of the Bethe-Salpeter equation corresponding to the vanishing energy-momentum of the 
electron-positron system is obtained. 
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I. INTRODUCTION 



It is well known that the structure of atoms (positronium included) is drastically modified 
by a magnetic field B if the field strength i? = IBI exceeds the characteristic atomic value 
B. ^ ,„Vc/.» . 2.35 X 10= G II ff. a .t.ong „,ag.et,c «e.d » B, tKe ™.a. 
perturbative treatment of the magnetic effects (such as Zeeman splitting of atomic energy 
levels) is not applicable, and instead, the Coulomb forces act as a perturbation to the 
magnetic forces. For positronium in such a field the characteristic size of the electron and 
positron across B is the Larmour radius 

Lb = (eB)-'/' = a,{BjBY'^ (1) 

and decreases with increase of the field strength, where Oq is the Bohr radius, Oq = (am)~^, 
a = 1/137. Henceforth, we set /?- = c = 1 and refer to the Heaviside-Lorentz system of units, 
where the fine structure constant a = e^/47r. 

The properties of positronium in a strong magnetic field {B ^ B^) are interesting for 
astrophysics because such fields are observed now for several kinds of astronomical compact 
objects (pulsars, powerful X-ray sources, soft gamma-ray repeaters, etc.). Besides, some of 
these objects are the sources of electron-positron pairs produced in their vicinities by various 
mechanisms jsf. At least a part of these pairs may be bound. For instance, at the surface 
of radio pulsars identified with rotation-powered neutron stars the field strength is in the 
range from ~ 10^ G to ~ 10^^ G j^j. A common point of all available^models of pulsars is 
that electron-positron pairs dominate in the magnetosphere plasma .5.]. These are formed 
by the single-photon production process in a strong magnetic field, 'j + B—>-e~^ + e~ + B. If 
the field strength is higher that ~ 4 x 10^^ G the pairs created are mainly bound 6]. Much 
more intense magnetic fields have been conjectured to be involved in several astrophysical 
phenomena. For instance, superconductive cosmic strings, if they exist, may have magnetic 
fields up to ~ 10"^^ — 10'^'^ G in their vicinities j?]. Electron-positron pairs may be produced 



near such strings 

In magnetic fields larger than i?a, the Coulomb force becomes more effective in binding 
the positronium because the charged constituents are confined to the lowest Landau level and 
hence to a narrow region stretching along the magnetic field (L^ ^ oq). Notwithstanding 
this effect, the binding energy of positronium AE is still very small in comparison with 
the rest mass, AE <^ 2m, even for the fields larger than Schwinger's critical value Bq = 



m^/e ~ 4.4 x 10^^ G, , i.e., the positronium remains an internally nonrelativistic system. 
The binding energy of the ground state, as calculated nonrelativistically, 

ma /. B ^ 



increases with increase of B, and the relativistic effects, for extremely huge fields, should 
be expected to become essential. The unrestricted growth of the binding energy (j2)) with 
the magnetic field is a manifestation of the fact that the Coulomb attraction force becomes 
supercritical in the one-dimensional Schrodinger equation, to which the nonrelativistic prob- 
lem is reduced in the high-field limit jl|, and the falling-to-the-center phenomenon occurs in 
the limit B = oo. 

Relativistic properties of positronium in a strong magnetic field were studied basing on 
the Bethe-Salpeter equation 0, [l^. The nontrivial energy dependence upon the transversal 



id, bj. Although the 



(pseudo)momentum component of the center-of-mass was found in 
Bethe-Salpeter equation is fully relativistic, it was used within the customary "equal-time" 
approximation that disregards the retardation effects, so that the relative motion of the 
electron and positron is treated in a nonrelativistic way. In this way the behavior Q is 

Q n 

reproduced for the ground state [21] - [lH . A completely relativistic solution for positronium 
in a strong magnetic field remains unknown. In this paper we study the positronium in 
an asymptotically strong magnetic field with not only the center-of-mass motion considered 
relativistically, but also the relative motion of its constituents. We point the ultimate value 
of the magnetic field guaranteeing such deepening of the positronium energy level that is 
sufficient to compensate for the whole rest mass 2m of it. 

To this end, in Section II we derive the fully relativistic - in two-dimensional Minkowsky 
space - form that the differential Bethe-Salpeter equation in the ladder approximation takes 
for the positronium when the magnetic field tends to infinity. This equation is efficient 
already for B ^ B^,. We also include a moderate external electric field parallel to B into 
this equation. In Section III the ultra-relativistic solution of maximum symmetry is found to 
the equation derived in Section II corresponding to the vanishing total energy-momentum 
of the positronium. The falling-to-the-center phenomenon |l^ characteristic of the two- 
dimensional equation of Section II for every positive value of the fine structure constant Q| 
is exploited for establishing the possibility that the zero energy point may belong to the 
spectrum, provided the magnetic field is sufficiently large. The origin of the falling to the 



center is in the ultraviolet singularity of the photon propagator. The effects of the mass 
radiative corrections and of the vacuum polarization are also considered. In concluding 
Section IV the results are summarized. 



II. TWO-DIMENSIONAL BETHE-SALPETER EQUATION FOR POSITRON- 
lUM IN AN ASYMPTOTICALLY STRONG MAGNETIC FIELD 

The view [3] that charged particles in a strong constant magnetic field are confined to 
the lowest Landau level and behave effectively as if they possess only one spacial degree 
of freedom - the one along the magnetic field - is widely accepted. Moreover, a conjecture 
exists jl^ that the Feynman rules in the high magnetic field limit may be directly served 
by two-dimensional (one space + one time) form of electron propagators. As applied to the 
^ 3ethe-Salpeter equation, the dimensional reduction in high magnetic field was considered in 
In these references the well-known simultaneous approximation to the Bethe-Salpeter 
equation taken in the integral form was exploited, appropriate for nonrelativistic treatment 
of the relative motion of the two charged particles. Once we shall in the next Section be 
interested in the ultrarelativistic regime, we reject from using this approximation, and find 
it convenient to deal only with the differential form of the Bethe-Salpeter equation. 

The electron-positron bound state is described by the Bethe-Salpeter amplitude (wave 
function) ')(x,f3ix:^, x^) subject to the fully relativistic equation (e.g., [lS|), which in the ladder 
approximation in a magnetic field may be written as 

= -iSnaD^^ix^ - xP)[7.]a/3[7,]m. X/3.(a;^ x^). (3) 

Here x*^, are the electron and positron 4-coordinates, Z)*-' (x'^— x^) is the photon propagator, 
and we have explicitly written the spinor indices A, (3, fi, u = 1,2,3,4. The metrics in the 
Minkowsky space is diag gij = {1, —1, —1, —1), i, j = 0,1,2,3. The derivatives 

d = 9^7, = + 5'7fc = 70^ + ^^^a^^' ^ = ^' 2' ^' 

act on x^ or x^ as indicated by the superscripts, and A = Aq'Jq — A^^jk- 

We consider the ladder approximation with the photon propagator taken in the Feynman 
gauge. With other gauges this approximation corresponds to summation of diagrams other 



than the ladder ones in agreement with the well-known fact that the ladder approximation 
is not gauge-invariant. 

We refer to, if needed, the so called spinor representation of the Dirac 7-matrices in the 
block form 





^01^ 




f -aA 


7o = 1 










I lOi 







(5) 



(Tfc are the Pauli matrices: 




102 




0-3 




(6) 



m is the electron mass, e the absolute value of its charge, e = 2^/twi. The vector potential of 
the constant and homogeneous magnetic field 5, directed along the axis 3 (-B3 = B, Bi^2 = 
0), is chosen in the asymmetric gauge 



AAx) 



Ao,2,3{x) = 0. 



(7) 



With this choice, the translational invariance along the directions 0,1,3 holds. 
Solutions to Eq. (jH)) may be represented in the form 



xP) = r]{xl - xl xl - xl x^^, ^2) exp{-[Po(x^ + x^) - P,{xl + x^)]}, (8) 



where Pq.s are the center-of-mass 4-momentum components of the longitudinal motion, that 
express the translational invariance along the longitudinal directions (0,3). Denoting the 



differences Xn 



Xn 



t, X' 



Xn 



from Eqs. ^ and (jH)) we obtain 



i9|| m + i9\ — 671^1 



A/3 



Pii 

-id\\ — - — m + i9^ + 671^1 (xg) 



X [r]{t,z,xY)]pu = -i87raAj(t,^,x^i^2 -<2) Na,3 bJ]^..[rj{t, z,xY)]pu , 
where x± = (xi, X2), d± = 71^ + 72^, ^ = ^70 + £73, and = Po7o - PsTs- 



(9) 



A. Fourier-Ritus Expansion in eigenfunctions of the transversal motion 

Expand the dependence of solution of Eq. Q on the transversal degrees of freedom into 
the series over the (complete set of) Ritus [l9| matrix eigenfunctions Eh{x2) 

h'^hP 
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Here Tjh°hp{^iZ) denote unknown functions that depend on the differences of the longi- 
tudinal variables, while the Ritus matrix functions e^^'^^^ Eh{x2) depend on the individual 
coordinates Xi \ transversal to the field. The Ritus matrix functions and the unknown func- 
tions Tfh'^hpit, z) are labelled by two pairs h^, of quantum numbers h = {k,pi, ), each pair 
relating to one out of the two particles in a magnetic field. The Landau quantum number k 
runs all nonnegative integers, k = 0,1, 2, 3..., while pi is the particle momentum component 
along the transversal axis 1. Recall that the potential A^{x) ((7j) does not depend on xi, so 
that pi does conserve. This quantum number is connected with the orbit center coordinate 
£2 along the axis 2 [12], pi = —X2eB. 

The matrix functions E'f^^ {X2) for transverse motion in the magnetic field ((Tj), re- 
lating in (jl(J|) to electrons (e) and positrons (p), are 4x4 matrices, formed, in the spinor 
representation, by four eigen-bispinors of the operator {—id± ± eA)^ 

i-id^ ± eA)l,e'P^^^[El'^X2)]i^-^^ = -2eBke'^^^^[El'^X2)]l^'^\ (11) 

placed, as columns, side by side Here the upper and lower signs relate to electron and 
positron, respectively, while a = ±1 and 7 = ±1 are eigenvalues of the operators 



i75 = , (12) 





diagonal in the spinor representation, to which the same 4-spinors are eigen-bispinors 

- i75i?r'^^ = lEt'^ , S34-^) = a^--) . (13) 

The couple of indices A = (a, 7) is united into one index A in the expansion (fTUj). A = 1, 2, 3, 4 
according to the convention: (+1,-1) = 1, (-1,-1) = 2, (+1,+1) = 3, (-1,+1) = 4. 
With this convention, the set of 4-spinors [Eh{x2\^i^''^^ = Eh{x2)^ can be dealt with as a 4 x 4 
matrix, the united index A spanning a matrix space, where the usual algebra of 7 -matrices 
may act. Correspondingly, in (fTUI) the unknown function [ffh^hpit, z)\\cxp is a matrix in the 
same space, and contracts with the Ritus matrix function. 

Following the matrix functions in expansion (jlUj) can be written in the block form 
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as diagonal matrices 



a^'P(/i; 0:1,2) 

a'='P(/i;xi,2) 



a°'n/i;xi,2)= "-'^ ' ' ^ I . (14) 

\ al\{h]Xi^2) 

Here a^'P(/i; xi,2) are eigenfunctions of the two (for each sign ±) operators 
[((— i(9_L)A ± eAxf =F creB] (we denote (i9±)a = d/dxx, A = 1, 2), labelled by the two values 
a = l, -1, 

[{{-id^)x ± eA^f T ae5]a:'P(/i; ^1,2) = 2eBka''f{h- x,,^), (15) 
namely, (we omit the subscript " 1" by pi in what follows) 

aTih; xi,2) = e'P^^U,^±^ [V^ (^0:2 ±^)], k = 0, 1, 2..., (16) 

with 

f/„(0 = exp |-^| (2"n!v/^)-i/7„(0 (17) 

being the normalized Hermite functions (iJ„(^) are the Hermite polynomials). Eqs. (fTH|) are 
the same as (jllj) due to the relation 

{A T eAf = -[(id^)x T eAxf ± e^Sg (18) 

and to Eq. (fT!^ . Simultaneously, the matrix functions (|T^ are eigenfunctions to the oper- 
ator — i(9i that commutes with E3 and 75 (fT^ . and with (i(9_L =F eA)^^,. The corresponding 
eigenvalue does not, however, appear in the r.-h. side of (fT3j) due to the well-known 
degeneracy of electron spectrum in a constant magnetic field. 
The orthonormality relation for the Hermite functions 

/oo 
Un{i)Un'{i)^i = 5nn'. (19) 
-oo 

implies the orthogonality of the Ritus matrix eigenfunctions in the form 

j El{x2)lEh'{x2f^dx2 = 6,k'Sxy. (20) 
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As a matter of fact, the matrix functions Eh{x2) are real, and we henceforth omit the complex 
conjugation sign 

The matrix functions e^^^^ E'^^^^ {X2) (jl4p commute with the longitudinal part 
± i(9|| — P\\/2 — m of the Dirac operator in Q, owing to the commutativity property 

Ii^^h(x2),7o,3]- = 0, (21) 



and are 19|, in a sense, matrix eigenf unctions of the transversal part of the Dirac operator 



(not only of its square (fTT|)) 

{id^ T ei)e'f^i£;^'P(x2) = ±^/2^ e'^^^ E^^ {x2)ii. (22) 

The Landau quantum number k appears here as a "universal eigenvalue" thanks to the 
mechanism, easy to trace, according to which the differential operator in the left-hand side 
of Eq. (j22|) acts as a lowering or rising operator on the functions (jl7j) . whereas the matrix 
o"2, involved in 72, interchanges the places the functions Uk^ Uk~i occupy in the columns. 
Contrary to relations, which explicitly include the variable a, whose value forms the number 
of the corresponding column, relations (fTT|) . (j22I); (PT|) . and the first relation in p3|) are 
covariant with respect to passing to other representation of 7-matrices, where the matrix 
Eh{x2) niay become non-diagonal. 



B. Equation for the Fourier-Ritus transform of the Bethe-Salpeter amplitude 



Now we are in a position to use expansion (fTUj) in Eq. Q. We left multiply it by 
(27r)-2e5 e-'P"^i£;|c(xi)e-'P''^?£;^p(x^), then integrate over d^x^a d^a;?,2- After using ^ 
and (pT|l . and exploiting the orthonormality relation piU) for the summation over the quan- 
tum numbers h^'^ = (A;'^'P,p^''^), the following expression is obtained for the left-hand side of 
the Fourier- Ritus-transformed Eq. Q: 



iSii 



Pii 



-7iV2e5F 



m — 7i 



J AA° 



Pll 



m + 7iV2e5fcP 



(23) 



We omitted the bars over the quantum numbers. 
Taking the expression 



E'lj (t, z, x-j^^ 



X 



9ij 
i47r2 



[X-, 



X 



P\2 



[Xo 



X 



p\2l-l 



(24) 



for the photon propagator in the Feynman gauge, we may then write the right-hand side of 
Ritus-transformed Eq. Q as 

a 



27r3 



^ + {xl - xiy + {xl - - t^' ' ^^^^ 
Integrating exphcitly the exponentials in ()25|) over the variable X = {xl + Xi)/2, we obtain 
the following expression: 

^ [dp dPi 6(P, - Pi) ^ g., I [El.{xlHEl{xl)U. [E'j^{xl)^A{xl)U. 

exp(ix(» — n))dx „ , „ -, r> 
^ ' z"^ + x"^ + {X2 — — t 

where the new integration variables x = xl — x^, Pi = + p^, p = {p'^ — p^)/2 and the 

new definitions Pi =p'^ + p^, p = {jf — P^)/2 have been introduced. The pairs of quantum 

numbers in ()2fj|l are 

^ ^pp^ Pi ^ Pi ^27) 

Hence the arguments of the functions (fTB|) in (^^1) are: 

^ (-^ - ^) ■ ^ (^^ - ^) ■ ^ ^) ■ ^ - ^) ■ 

successively as the functions Eh{xi^2) appear in (j^Uj) from left to right. After fulfilling the 
integration over dPi with the use of the 5-function, introduce the new integration variable 
q = p — p instead of p, and the integration variables x^ = x^ + {Pi + 2p)/2eB, X2 = 
X2 — {Pi — 2p)/2eB instead of x% and x^. Then may be written as 



^ / dg ^ g,, l[El.{xl)^,El.{x'2)U^[E2^{xl)r^^^^^ 

k^kP 



AoAP 



^ . exp(— ixg) dx eB dxj dx^ ^^g-j 



^2 + ^2+ {xl -xl- ^) - t2 
Now the pairs of quantum numbers in (j29p are 

^e,p ^ ^^,p^ Pi ^ ^ ^ ± g ± p). (30) 

Hence the arguments of the functions ()16p in ()29p from left to right are 

yipr2, yip(^r2 + ^), ^{xi--^y (v^xp). (31) 



C. Adiabatic approximation 



Now we aim at passing to the large magnetic field regime in the Bethe-Salpeter equation, 
with (j23p as the left-hand side and (j29j) as the right-hand side. Define the dimensionless 
integration variables w = x^feB, q' = q/^feB, ^'^'P = li%f\/eB in function (|29|) . Then it 
takes the form 



A<=AP 



^ I exp(-iwgOdwd^M^P ^^^^ 



The pairs of quantum numbers in ()32|1 are 



h'" = {k ^ ± P), /^^'P = Y ± ^VeP ± p). (33) 
The arguments of the functions (fTBj) in (jH^ from left to right are 

r, e + q, e-q, e- m 

When considering the large field behavior we admit for completeness that the difference 
between the centers of orbits along the axis 2 X2 — x^ = —j^ may be kept finite, in other 
words that the transversal momentum Pi grows linearly with the field. We shall see that 
that big transversal momenta do not contradict dimensional compactification, but produce 
an extra regularization of the light-cone singularity. 

In the region, where the 2-interval {z"^ — t^y/"^ essentially exceeds the Larmour radius 
Lb = (eP)-^/^ 

^2 - ^2 ^ L| , (35) 

one may neglect the dependence on the integration variables w and later on ^^'"^ in the 
denominator. Integration over w produces 2Tr6{q'), which annihilates the dependence on q' 
in the arguments (jHT|) of the Hermite functions, and they all equalize. 

Let us depict this mechanism in more detail. Fulfill explicitly the integration over dw in 



exp(— iwg') dw y/eBn 



X 



9{q') exp [-gVe5(^' -^') + ^'] + ^("^O exp [ q ^ eB{z'' - t^) + ] , (36) 
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where 



A'=[c-e-4k-^') (37) 



and 9{q') is the step function, 



1 when q' > 0, 

^W) = I \ when q' = 0, (38) 
when q' < 0. 

Due to the decreasing exponential in (fTTjl the variables ^'^'P do not exceed unity in the 
order of magnitude and can be neglected as compared to in (jHTj). Unless q' is large it 
may be neglected as compared to the same term in (jHTjl . too. Then = and after 

()36|) is substituted in ()32j) and integrated over dg' the contribution comes only from the 

1 

integration within the shrinking region \q'\ < {eB[z'^ — + Then q' can be also 

neglected in the arguments (jHll). If, contrary to the previous assumption, we admit that 
\q'\ is of the order of ■^^= ~ V eB we see that the exponentials in (jHBj) fast decrease with 
the growth of the magnetic field as exp [— ei?(2;^ — t^)], and therefore such values of \q'\ do 
not contribute to the integration. If we admit, last, that \q'\ ^ l^^l^ we find that the 



contribution exp [—\q' \ ^ye^iz^~—W)~+J/f)^ ] from the integration over such values is still 
smaller. Thus, we have justified the possibility to omit the dependence on q' in (jTTjl and in 
(I34p . and also on ^'^'P in (j37j) . Now we can perform the integration over dg' to obtain the 
following expression for 



2 I Pf ,2 XI 



It remains yet to argue that the limit (j39j) is valid also when the term ^ is not kept. In 
this case we no longer can disregard q' inside when q' is less than or of the order of unity. 
But we can disregard as compared with eB{z'^ — t^) to make sure that the integration 
over dg' is restricted to the region close to zero \q\' < [eB{z^ — and hence set g' = 
in (jMj) . The contribution of large q' is small as before. 

The integration over ^'''P of the terms with i = 0, 3 in (j39p yields the Kroneker deltas 
^kPkp orthonormality (fT^ of the Hermite functions thanks to the commu- 

tativity (PT|) of the Ritus matrix functions with 70 and 73. On the contrary, 71, 72 
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do not commute with ()14|) . This imphes the appearance of terms, non-diagonal in Landau 
quantum numbers, hke Sf^c-j^^i and 6|^p-gl>^-^^, in (jH^ . proportional to (z = 1,2): 

-A 



E 



.° ^'"'^ I [VhMt,z)l.x.. (40) 
A- n 



Here expressed in terms of ^ through the chain of the changes of variables made 

above starting from (j2Sl), so that all the arguments of the Hermite functions have become 
equal to C,. Besides, 



h^'P = (r'P,p"'P), = (F'^p^'P), / +pP = Pi. 



(41) 



'^L = / aih,X2) (Ji a{h,X2)d^, i = l,2. 



(42) 



kk 



a'_i{h, X2)a'_^_i{h, X2) 



a'j^^{h,X2)a'_^{k,X2) 




& 



(43) 



aP) = i 

kk 







-a+i(/i,X2)a_i(A;,X2) \ . / -5- 

I = 1 





k, k-1 

St u,-, 



(44) 



a'_^{h,X2)a'_^i{h,X2) u y 

The prime over a indicates that the exponential exp(ipxi) is dropped from the definitions ()14|) 
and (|T6|l . The non-diagonal Kronecker deltas appeared, because a'_^i{h,X2) are multiplied 
by a'^i{h,X2) under the action of the cri^2-blocks in 71^2 dS}- In the final form, the matrices 
in (jini) are 



-A 



kk 



Ai 

kk 



2 (7i(±%fe-i + h,k+i) + h2{±Sk,k-i - h,k+i)) 



(45) 



with the upper sign relating to i = 1 and the lower one to i = 2. Now Eq. ^ acquires the 
following form: 



idn — m — ^iV2eBk'^ 

2 



AA= 



-i(9|| — m + 7i V2eBkP 



A<=AP 



^AP 



p2 

^ + (eB)2 V=0,3 



t2 



XI ^u[7i]AA47*]^aP [?7/^o;,p(t, z)]xoXP - Yl ^kLl, fcp±i • (46) 



i=l,2 
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The bars over quantum numbers are omitted. This equation is degenerate with respect to 
the difference of the electron and positron momentum components p = {p^ —p^)/2 across the 
magnetic field, but does depend on its transversal center-of-mass momentum Pi = {jf + p^). 
This dependence is present, however, only for sufficiently large transverse momenta Pi. 

At the present step of adiabatic approximation we have come, for high magnetic field, to 
the chain of Eqs. ()46|) . in which the unknown function for a given pair of Landau quantum 
numbers fc*^, kP is tangled with the same function with the Landau quantum numbers both 
shifted by ±1 (in contrast to the general case of a moderate magnetic field, where these 
numbers may be shifted by all positive and negative integers). To be more precise, the chain 
consists of two mutually disentangled sub-chains. The first one includes all functions with 
the Landau quantum numbers /c*^, kP both even or both odd, and the second includes their 
even-odd and odd-even combinations. We discuss the first sub-chain since it contains the 
lowest function with k!^ = kP = 0. We argue now that there exists a solution to the first 
sub-chain of Eqs. (j46|l . for which all T/^o po. fcp,pp(t, z) disappear if at least one of the quantum 
numbers /c*^, k^ is different from zero. Indeed, for k'^ = k^ = Q Eqs. then reduces to the 
closed set 

= -^—Jq 5Z 9iibi\\Al'i\t^\-^ iVo^pi-Ap^lit, z)]x^xp . Pi + = ^1- (47) 

In writing it we have returned to the initial designation of the electron and positron trans- 
verse momenta p1'^. Denote for simplicity ?7fc°fcp =''7fc^Pl;fcp,pf -^)- consider Eqs. fHUj) 
with k'^ = k^ = 1 for Tju we shall have a nonzero contribution in the right-hand side, pro- 
portional to 7^00 coming from T^c_i^p_i, since the other contributions Tj ii,Tj 22, 'ij 20, IJ02 are 
vanishing according to the assumption. As the left-hand side of Eq. (j46j) now contains a 
term, infinitely growing with the magnetic field B, it can be only satisfied with the function 
7^11, infinitely diminishing with B in the domain as 

[Vn]xf^ = ~ 2e^ 2 % 2 5Z 9ii[lili]xAlili],^x^ [T/ooIa-ap . (48) 

in accord with the assumption made. Thus, the assumption that all Bethe-Salpeter ampli- 
tudes with nonzero Landau quantum numbers are zero in the large-field case is consistent. 
We state that a solution to the closed set (jTrf) for TjQpo.Qpp(t, z) with all the other components 
equal to zero is a solution to the whole chain pUjl. 

13 



The derivation given in this Subsection reahzes formally the known heuristic argument 
that, for high magnetic field, the spacing between Landau levels is very large and hence the 
particles taken in the lowest Landau state remain in it. Effectively, only the longitudinal 
degree of freedom survives for large B, the space-time reduction taking place. Eq. (j47j) is 
a fully relativistic two-dimensional set of equations with two space-time arguments t and 



z and two gamma-matrices 70 and 73 involved. Since, unlike the previous works 
11 1 , neither the famous equal-time Ansatz for the Bethe-Salpeter amplitude nor any 
other assumption concerning the non-relativistic character of the internal motion inside the 
positronium atom was made, the equation derived is valid for arbitrary strong binding. It 
will be analyzed for the extreme relativistic case in the next Section. 

The two-dimensional equation (j47j) is valid in the space-like domain (j35|) . It is mean- 
ingful provided that its solution is concentrated in this domain. In non-relativistic or semi- 
relativistic consideration it is often accepted that the wave function is concentrated within 
the Bohr radius = (am)~^ ^ 0.5 x 10^*^ cm. It is then estimated that the corresponding 
analog of asymptotic equation (jTfjl holds true when qq ^ Lb, i.e. for the magnetic fields 
larger than = oP'm^ je ~ 2.35 x 10^ G. This estimate, however, cannot be universal and 
may be applicable at the most to the magnetic fields close to the lower bound where the value 
of the Bohr radius can be borrowed from the theory without the magnetic field. Generally, 
the question, where the wave function is concentrated, should be answered a posteriori by 
inspecting a solution to Eq. (j47j) . Therefore, one can state, how large the fields should be 
in order that the asymptotic equation (jUj) might be trusted, no sooner that its solution is 
investigated. We shall return to this point when we deal with the ultra-relativistic situation. 

Remind that the transverse total momentum component of the positronium system is 
connected with the separation between the centers of orbits of the electron and positron 
Pi/{eB) = X2 — 5^ in the transversal plane, so that the "potential" factor in Eq. ()47p may 
be expressed in the following interesting form 

(49) 



{cf the corresponding form of the Coulomb potential in the semi-relativistic treatment of 
the Bethe-Salpeter equation in [l^ 11 1- the difference between the potentials in Q, [ll| 
lies within the accuracy of the adiabatic approximation). The appearance of in the 
potential determines the energy spectrum dependence upon the momentum of motion of the 
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two-particle system across the magnetic field like in 2ll |. 

We shall need Eq. ()47|) in a more convenient form. First, transcribe it as 

- ^ - m) 77o,p-o,pf (i^, z) (-i5|| - ^ - m)T 

" 2^ ^ ^o,p-o,p?(i, ^) 7i • (50) 

Here the superscript T denotes the transposition. With the help of the relation '~^J = 
—C^^'-fiC, with C being the charge conjugation matrix, = 1, and the anti-commutation 
relation [74,75]+ = 0, 7I = —1, we may write for a new Bethe-Salpeter amplitude Q(t,z), 



defined as 



the equation 



0(^>^)= r/o,pC;0,p?(t,^)C75, (51) 



~^ P\\ ^ P\\ 

iA\ - y -"^)0(^'^)(-i^ii - Y - H 



= ^^^S galiQit, z)'yi . (52) 

^ + (eB)2 i=0,3 

The unknown function here is a 4x4 matrix, which contains as a matter of fact only 
four independent components. In order to correspondingly reduce the number of equations 
in the set ()52|) . one should note that the 7- matrix algebra in two-dimensional space-time 
should have only four basic elements. In accordance with this fact, only the matrices 70,3 
are involved in ()52|) . Together with the matrix 7073 and the unit matrix I they form the 

basis, since 70,3 ■ 7o73 = 73,o, 7o = "71 = (7073)^ = 1, [70,73]+ = [7o,3,7o73]+ = 0. Using 
this algebra and the general representation for the solution 

e = a/ + 670 + C73 + (i7o73, (53) 

one readily obtains a closed set of four first-order differential equations for the four functions 
a, b, c,d of t and z. The same set will be obtained, if one replaces in Eqs. (fS^ and (fS^ the 
4x4 matrices by the Pauli matrices (jHI), subject to the same algebraic relations, according, 
for instance, to the rule: 70 ^ cr3, 73 ^ io"2, 7073 =^ (Ji- Then Eq. f|7fjl becomes a matrix 
equation 

^ ^ p p. i i P P 

{idtcr^ + 9^0-2 - ^0-3 + ^10-2 - m){}(t, z){-idt(T3 - 9^(72 - —cr^ + ^1(72 - m) 



^""^ [a^§{t, z)a-i + (T2^9(t, z)a2] (54) 



z^ H — - — 



(eB) 
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for a 2x2 matrix ^9, 

d = al + ha^ + 1072 + dai . 
Here / is the 2x2 unit matrix, and functions a, 6, c, ci are the same as in ()53|) . 



(55) 



D. Including an external electric field 



Let us generahze the two-dimensional Bethe-Salpeter equation obtained in the presence 
of a strong magnetic field by including an external electric field, parallel to it, that is not 
supposed to be strong, E <^ B. To this end we supplement the potential ((Zj) in Eq. (jHl) by 
two more nonzero components 



^0(2:0,2:3), ^3(2:0, 2:3)) 7^ 0, 



(56) 



that carry the electric field - not necessarily constant - directed along the axis 3. We shall use 
the collective notations A\\ = {Aq, A3), x\\ = (xq, xs), Sy'^ = dQ^-yo — dl'^'js, Ay = ylo7o— ^4373. 
We shall not exploit now a representation like (jHl), but deal directly with the Bethe-Salpeter 
amplitude )({x^,x^) as a function of the electron and positron coordinates, and with its 

e P 

"ll'^ll 

(fTn|) . In place of Eq. one should write 



Fourier-Ritus transform ')(h'=hp{x?,,x?:) connected with ')(^{xf,,xfr,x\,x'^) in the same way as 



idn 



ei||(x|j) -m + idl- e-fiAiix^) 



J A/3 



idn 



eij|(x?) ~m + idl + 671^1 (x^) 



X [x(xJ'P, xr)]pu = -i8vr«A,(t, z, xl^ - x^^) [7.]/.. [x{^Y . ^7)1/3. • (57) 



Thanks to the commutativity ()2ip the rest of the procedure of the previous Subsection 
remains essentially the same, and we come, in place of (07j), to the following two-dimensional 
equation: 



\d\\ — eA\\ ixl 



- m 

J A/3 



id\\ +eA||(X|| 



m 



J 



Xo,p5;0,p?(2^|j) 3;n)J/3i/ 



P? 



9ii [lilxa [li] 



(58) 



(eB)2 «=0,3 

for a positronium atom in a strong magnetic field placed in a moderate electric field, parallel 
to the magnetic one. In order to apply this equation to a system of two different oppositely 
charged particles interacting with each other through the photon exchange and placed into 
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the combination of a strong magnetic and an electric field in the same direction, say a 
relativistic hydrogen atom, one should only distinguish the two masses in the first and 
second square brackets in the left-hand side. 



III. ULTRA-RELATIVISTIC REGIME IN A MAGNETIC FIELD 

In the ultra- relativistic limit, where the positronium mass is completely compensated by 
the mass defect, Pq = 0, for the positronium at rest along the direction of the magnetic field 
P3 = 0, the most general relativistic-covariant form of the solution (j53|) is 

e = /$ + 9||$2 + 7o73<^'3- (59) 

The point is that 7073 is invariant under the Lorentz rotations in the plane {t,z). Substi- 
tuting this into ()52|) with Pq = -P3 = we get a separate equation for the singlet component 
of §^ 

-□2 + m") <l>(t, z) = J- ' (60) 

and the set of equations 

4a7r"i$3(t,z) 



(□s + m") $3(1 



z 



El 

(eB) 

2 



(-□2 + m')dt^2 + 2mi9,<l>3 = 0, 
(-□2 + m^)d,^2 + 2mi9i<l>3 = (61) 

for the other two components. Here ^2 = —d'^/dt^ + d"^ /dz^ is the Laplace operator in two 
dimensions. Note the "tachyonic" sign in front of it in the first equation ()61|). 

Let us differentiate the second equation in (jUT|) over z and the third one over t and 
subtract the results from each other. In this way we get that □2^3 = 0. This, however, 
contradicts the first equation in (|F)T|l if $3 7^ 0. Therefore, only $3 = is possible. Then, 
the two second equations in (jHTjl are satisfied, provided that (—02 + m^)$2 = 0. We shall 
concentrate in Eq. ()60|) in what follows. 

The longitudinal momentum along Xi, or the distance between the orbit centers along X2, 
plays the role of the effective photon mass and a singular potential regularizator in Eq. ()6Up . 
The lowest state corresponds to the zero value of the transverse total momentum Pi = 0. In 
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this case Eq. ()60p for the Ritus transform of the Bethe-Salpeter amphtude finally becomes 

(-□2 + m^)$(t,^) = -^-^. (62) 

We consider now the consequences of the fall-down onto the center phenomenon present in 
Eq. ()62|) . formally valid for an infinite magnetic field, and the alterations introduced by its 
finiteness. 



A. Fall-down onto the center in the Bethe-Salpeter amplitude for strong magnetic 
field 

In the most symmetrical case, when the wave function $(a;) = does not depend 
on the hyperbolic angle in the space-like region of the two-dimensional Minkowsky space, 



t = s sinh 0, 2; = s cosh 4>, s = \/ — f^, Eq. (jU^ becomes the Bessel differential equation 

d2$ ld$ 2^ 4a ^ 
- -— - -— + m2$ = — $. 63 

ds^ s as TTs"^ 

It follows from the derivation procedure in the previous Section II that this equation is valid 
within the interval 

< So < s < 00, (64) 



where the lower bound sq depends on the external magnetic field - it should be larger than 
the Larmour radius Lb = (eB)^^^^ and tend to zero together with it, as the magnetic field 
tends to infinity. The stronger the field, the ampler the interval of validity, the closer to the 
origin s = the interval of validity of this equation extends. If the magnetic field is not 
sufficiently strong, the lower bound sq falls beyond the region where the solution is mostly 
concentrated and the limiting form of the Bethe-Salpeter equation becomes noneffective, 
since it only relates to the asymptotic (large s) region, while the rest of the s-axis is served 
by more complicated initial Bethe-Salpeter equation, not reducible to the two-dimensional 
form there. This is how the strength of the magnetic field participates - note, that the 
coefficients of Eq. (jHSl) do not contain it. 

Solutions of (jU^ behave near the singular point s = like s'^, where 



a = ±2j-^. (65) 
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The fall-down onto the center 
the genuine value a - 
field where = vr/ 



12| occurs, if a > = 0, i.e., for arbitrary small attraction, 
-/137 included. This differs crucially from the case of zero magnetic 
isl ]. This difference is a purely geometrical consequence jl^ of the 
dimensional reduction of the Minkowsky space from (1,3) to (1,1). 

In discussing the physical consequences of the falling to the center we appeal to the 



approach recently developed by one of the present authors as applied to the Schrodinger 



equal 



field 



ion with singular potential 



22| and to the Dirac equation in supercritical Coulomb 



23|. Within this approach the singular center looks like a black hole. The solutions of 



the differential equation that oscillate near the singularity point are treated as free particles 
emitted and absorbed by the singularity. This treatment becomes natural after the differen- 
tial equation is written as the generalized eigenvalue problem with respect to the coupling 
constant. Its solutions make a (rigged) Hilbert space and are subject to orthonormality 
relations with a singular measure. This singularity makes it possible for the oscillating 
solutions to be normalized to (5-functions, as free particle wave-functions should be. The 
nontrivial, singular measure that appears in the definition of the scalar product of quantum 
states in the Hilbert space of quantum mechanics introduces the geometry of a black hole of 
non-gravitational origin and the idea of horizon. The deviation from the standard quantum 
theory manifests itself in this approach only when particles are so close to one another that 
the mutual Coulomb field they are subjected to falls beyond the range, where the standard 
theory may be referred to as firmly established [2^. 

Following this theory we shall be using sq as the lower edge of the normalization box 
I23I. For doing this it is necessary that sq be much smaller than the electron Compton 
length. So <^ ~ 3.9 x 10^^^ cm, the only dimensional parameter in Eq. ()63|). In this case 
the asymptotic regime of small distances is achieved and nothing in the region s < sq beyond 
the normalization volume - where the two-dimensional equations (jlTj) . (jSH), (p?)). (IHH) and 
hence ()63|1 are not valid and the space-time for charged particles remains four-dimensional 
- may affect the problem, because this is left behind the event horizon. 

In alternative to this, we might treat sq as the cut-off parameter. In this case we have 
had to extend Eq. continuously to the region < s < Sq, simultaneously replacing 

the singularity s'"^ in it by a model function of s, nonsingular in the origin, say, a constant 
Sq^. In this approach the results are dependent on the choice of the model function which 
is intended to substitute for the lack of a treatable equation in that region. Besides, the 
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limit So ^ does not exist. The latter fact implies that the approach should become invalid 
for sufficiently small sq, i.e., large B. We, nevertheless, shall also test the consequences of 
this approach later in this section to make sure that in our special problem the result is not 
affected any essentially. 



B. Ultimate magnetic field 

With the substitution = '^{s)/y/s Eq. (jUHj) acquires the standard form of a 

Schrodinger equation 

- + % ' ^(g) + ^'^(s) = . (66) 

Equation ()66|) is valid in the interval 

So < s < oo. So > Lb = (e5)"^/2 (67) 

Treating the applicability boundary sq of this equation as the lower edge of the normalization 
box, as discussed above, sq -C m"^, we impose the standing wave boundary condition 

^{so) = 0, (68) 

on the solution of (I^Uj) 



^(s) = JCuims), V = 12 V«A - 0.096 i (69) 
that decreases at infinity. It behaves near the singular point s = as 

.2) r{i + u) ~ (2) r{i-u)' ^''^^ 

Here the Euler F-functions appear. Starting with a certain small value of the argument 
ms, the McDonald function with imaginary index K,jj{ms) ()69|1 oscillates, as s — > 0, passing 
the zero value infinitely many times. Therefore, if Sq is sufficiently small the standing wave 
boundary condition (jUHj) .prescribed by the theory of Refs. 2^, 13' t)e definitely satisfied. 
Keeping to the genuine value of the coupling constant a = 1/137 (z/ = 0.096 i) one may ask: 
what is the largest possible value s™^^ of Sq, for which the boundary problem (j66|) . (j68j) can 
be solved? By demanding, in accord with the validity condition of Eqs. and 
that the value of s^^^ should exceed the Larmour radius: 

s— »(ei?)-V2 B»-^^^, (71) 
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one establishes, how large the magnetic field should be in order that the boundary problem 
might have a solution, in other words, that the point Pq = P = might belong to the 
spectrum of bound states of the Bethe-Salpeter equation in its initial form 0. 

One can use the asymptotic form of the McDonald function near zero to see that the 
boundary condition (j68p is satisfied provided that 



/mso\2- _ r(l + z^ 

\ 2 ) ~ T*(l-u) ^^^^ 



or 



z/ln^ = iargr(z/ + 1) -ivrn, n = 0, ±1, ±2, ... (73) 

Since is small we may exploit the approximation r(l + z/) ^ 1 — z/Ce, where Ce = 0.577 
is the Euler constant, to get 

%3 



We have expelled the non-positive integers n from here, since they would lead to the roots for 
mso of the order of or larger than unity in contradiction to the adopted condition sq <^ m~^. 
For such values eq. ()7()j) is not valid. It may be checked that there are no other zeros of 
McDonald function, besides (f7i|) . The maximum value for sq is provided by n = 1. We 
finally get 

.(^)^-i;?-c. 

or 

= -exp j-iy?- Ce j ^ 10--1. (75) 
m I 2 V a J m 

This is fourteen orders of magnitude smaller than the Compton length = 3.9 x 10~^^cm 
and makes about 10^^^ cm. Now, in accord with (|7T| . if the magnetic field exceeds the 
ultimate value of 

^2 ( 3/2 >, 

= — exp I + 2Ce [ ^ 1.6 X 10^8 Bo, (76) 



the positronium ground state with the center-of-mass 4-momentum equal to zero appears. 
Here Bq = m?/e ~ 1.22 x 10^^ Heaviside-Lorentz units is the Schwinger critical field, or 
Bq = w?c^/eh ~ 4.4 X 10^^ G. The value of B^i^ is ~ 10"'^ G that is a few orders of 

21 



magnitude smaller than the highest magnetic field in the vicinity of superconductive cosmic 
strings [ij. Excited positronium states may also reach the spectral point = 0, but this 
occurs for magnetic fields, tens orders of magnitude larger than (j76p - to be found in the 
same way from (j74j) with n = 2, 3... 

The ultra-relativistic state = has the internal structure of what was called a confined 
state in j^;!^' i-^- whose wave function behaves as a standing wave combination 

of free particles incoming from behind the lower edge of the normalization box and then 
totally reflected back to this edge. It decreases as exp(— ms) at large distances like the wave 
function of a bound state. The effective "Bohr radius", i.e. the value of s that provides 
the maximum to the wave function (|69|) makes Smax = 0.17m^^ (this fact is established by 
numerical analysis). This is certainly much less than the standard Bohr radius Oq = {am)~^. 
Taken at the level of 1/2 of its maximum value, the wave-function is concentrated within 
the limits 0.006 m^^ < s < 1.1 m~^. But the effective region occupied by the confined state 
is still much closer to s = 0. The point is that the probability density of the confined state is 
the wave function squared weighted with the measure s~^ds singular in the origin 2^, [2^ 
and is hence concentrated near the edge of the normalization box sq ~ 10^^^ cm, and not 
in the vicinity of the maximum of the wave function. The electric fields at such distances 
are about 10^^ Volt/cm. Certainly, there is no evidence that the standard quantum theory 
should be valid under such conditions. This remark gives the freedom of applying the theory 
presented in Refs. Q> Q- 

A relation like (f7F)|l between a Fermion mass and the magnetic field is present in [2^ . 
There, however, a different problem is studied and, correspondingly, a different meaning 
is attributed to that relation: it expresses the mass acquired dynamically by a primarily 
massless Fermion in terms of the magnetic field applied to it. The mass generation is 
described by the homogeneous Bethe-Salpeter equation, whose solution is understood 13,01 
as the wave function of the Goldstone boson corresponding to the spontaneous breaking of 
the chiral symmetry characteristic of the massless QED. It is claimed, moreover, that the 
resulting relation between the magnetic field and the acquired mass is independent of the 



choice of the gauge for the photon propagator. The equations of Ref. [2J| may well be read 
off, formally, as serving our problem of the compensation of the positronium rest mass by the 
mass defect in a magnetic field, too, and the resulting expression may be used for determining 
the corresponding magnetic field, provided that the electron mass m is substituted for the 
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acquired mass mdyn of 2J|. There is, however, an important discrepancy in numerical 



coefficients in the characteristic exponential between ()76p and the corresponding formula in 



2J: the latter contains exp{n^^^ / {2ay/'^} in place of exp{7r^/V«^^^ + 2Ce} in ^ and its 



direct use would lead to a more favorable estimate of the ultimate value of the magnetic 



field, 2.6 x lO^^^o, than Although the 

and falling to the center, acting here and in 



Dasic mechanisms, the dimensional reduction 
2J], are essentially the same, the procedures 



are very much different, and the origin of the discrepancy remains unclear. Later, in the 
authors revised their relation in favor of a different approximation. Supposedly, the revised 
relation may be also of use in the problem of ultimate magnetic field dealt with here. 

It is interesting to compare the value (I76| with the analogous value, obtained earlier 
by the present authors (see p. 393 of Ref. [lO|) by extrapolating the nonrelativistic result 
concerning the positronium binding energy in a magnetic field to extreme relativistic region: 

i?ultUoN-REL = ^^-P - lO'^'^O. (77) 

Such is the magnetic field that makes the binding energy of the lowest energy state equal to 
(-2m). (This is worth comparing with the magnetic field, estimated j2^ as a'^ exp{2/a)Bo, 
that makes the mass defect of the nonrelativistic hydrogen atom comparable with the elec- 
tron rest mass). We see that the relativistically enhanced attraction has resulted in a 
drastically lower value of the ultimate magnetic field. Note the difference in the character 
of the essential nonanalyticity with respect to the coupling constant: it is exp{'^^^/7l/^/a) in 
fl76l) and exp(2v^/Q;) in (1771) . Another effect of relativistic enhancement is that within the 

□ n 

semi-relativistic treatment of the Bethe-Salpeter equation |2l|-|lH, as well as within the one 
using the Schrodinger equation Q], only the lowest level could acquire unlimited negative 
energy with the growth of the magnetic field, whereas according to (|74|) in our fully rela- 
tivistic treatment all excited levels with n > 1 are subjected to the falling to the center and 
can reach in turn the point Py =0. 

Let us see now, how the result (fTBj) is altered if the cut-off procedure of Ref. jl^ is used. 
Consider Eq. (j66|l in the domain sq < s < oo, but replace it with another equation 

- - ^^Ms) + m'<fo{s) = (78) 

as Sq 

in the domain < s < sq. The singular potential is replaced by a constant near the origin 
in dZHl). Demand, in place of that ^o(O) = 0, [^[,(so)/^o(so)) = (^'(so)/^(so)]- Then, 
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the result ()76|) will be modified by the factor 

exp I == I , (79) 

I A/¥TIcot(^ + i)-lJ 

which may be taken at the value a = 0. Thus, the result ()76|) is only modified by a factor of 
exp(— 4/3) ^ 0.25. Generally, the estimate of the limiting magnetic field (fTHj) is practically 
nonsensitive to the way of cut-off, in other words to any solution of the initial equation 
inside the region < s < sq, where the magnetic field does not dominate over the mutual 
attraction force between the electron and positron. This fact takes place, because the term 
(tt^/^/y^) ^ 65, singular in a, is prevailing in ()76|). the details of the behavior of the wave 
function close to the origin s = being not essential against its background. 



C. Radiative corrections 

1. Vacuum polarization 

We should answer the question of whether the effects of vacuum polarization in a strong 
magnetic field may or may not screen the interaction between the electron and positron 
in such a way as to prevent the falling to the center in the positronium atom. It is clear 
aposteriory that no matter how strong the magnetic field is, the ultraviolet singularity dom- 
inates over its infiuence in the photon propagator, if the interval sufficiently close to the 
light cone is involved. Therefore, there is a competition between the magnetic field and this 
characteristic interval, which is in our problem the Larmour radius that itself depends on 
the magnetic field. We have to consider the outcome of this competition. 

To include the effect of the vacuum polarization we should use the photon propagator in 
a magnetic field, whose infiuence is realized via the vacuum polarization radiative correc- 
tions, instead of its free form ()24|) used above. The photon propagator in a constant and 
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lomogeneous magnetic field has the following approximation-independent structure 

D^Jix) = j exp{ikx)Di,ik) d% t,j = 0, 1,2,3, (80) 
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arbitrary, 



a = 1,2,3 



a=4 



(81) 



Here ft*-"^ and are four eigenvectors and four eigenvalues of the polarization operator lljj 

U/ hf = K,{k) hf\ (82) 



The eigenvectors are known in the final form: 



(1) 



{F''k),e - k,{kF''k), 6f = (FA;) 



l{3) 



(Fk), 



'(4) 



ki. 



(83) 



where F, F and are the external electromagnetic field tensor, its dual, and its tensor 
squared, respectively, contracted with the photon 4-momentum k. On the contrary, the 
eigenvalues Ki^2,3{k) are generally unknown - subject to approximate calculations - scalar 
functions of two Lorentz-invariant combinations of the momentum and the field, which in 
the special frame, where the external electromagnetic field is given by ((7j), are fcg — fc| and 
kf + k2 = k\. The eigenvalue K4 is equal to zero as a trivial consequence of the transversality 
Tl^kj = of the polarization operator. The eigenvectors (j83|l with a = 1, 2, 3 are 4-potentials 
of the three photon modes, while the dispersion laws of the corresponding electromagnetic 
eigenwaves are obtained by equalizing the denominators in ()8H1 with zero. In the special 
frame the eigenvectors (jH^ . up to normalizations, are 



/o\ 

ki 
k2 



+ k] 



fko\ 

ki 
k2 



l(2) 



fkA 





(,(3) 



/ \ 

k2 
-ki 

\ y 



(84) 



When calculated j^^], [2^ within the one-loop approximation of the Furry picture (i.e. 
using exact Dirac propagators in the external magnetic field without radiative corrections) 

nnnnn 

these eigenvalues have the following asymptotic behavior [30|, [28[, [29|, [31[, [32| (note the 
difference in the signs in front of fc^ due to a different metric convention used here) for large 



fields eB > m^, eB > \kl - k^ 



Ki(kQ k^,kj 



—ak"^ 
Sir 



1.21 



(85) 
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TiBo V '^m? B ) 4m2 - {kl - A;|)(l - r/2) 



— aA;^ f , B \ a 



(86) 



3. (>n|-c)-^[0.2Ui-1.21(.= -*,^)]. (87) 

C = 0.577 is the Euler constant. Eqs. (jHH|) and (jHTj) are accurate up to terms, decreasing 
with B hke {Bcr/B)\n{B/ Bcr) and faster. Eq. (jHEI) is accurate up to terms, logarithmically 
growing with B. In we took also the limit k']_ ^ (B / Bcr)iTL^ , which is not the case for 
K2, wherein the factor exp {—k\BQ/2m^B) is kept different from unity. Although the 
components ^1,2,3 contain the growing logarithms a \n{B/Bo) the latter are yet small for the 
values of the magnetic field of the order of i?uit ()76|). This is not the case for the linearly 
growing part of 

Let us inspect the contributions of the photon propagator (jHT|l into the equation that 
should appear in place of ()47|) . To match the diagonal form ()24j) corresponding to the 
Feynman gauge we fix the gauge arbitrariness by choosing 

D^{k) = + Ki{k)]-\ (88) 

In the isotropic case where no magnetic field is present all the three nontrivial eigenvalues 
are the same, Ka{k) = ^{k), a = 1,2,3. Then, with the choice (j88p in (j8H) the photon 
propagator in this limit becomes diagonal 

^''^^^ ^ ~k^ + K{k) ^ lu-yf ^ ~k^ + K{ky ^^^^ 

since the eigenvectors (jH^ or (jS^ make an orthogonal basis irrespective of whether the 
magnetic field is present or not. 

In spite of the presence [s^ of a term, linearly growing with the field in (j86|) . the 
component D2 does contribute in the limit of high fields into the right-hand side of an 
equation to replace (jUUj) . because the ultra-violet singularity at the distance of the Larmour 
radius from the light cone dominates. To see this note that the right-hand side of the analog 
of (|T7|) should get the contribution from D2: 

d^fc . (90) 



1 f [hlo - tolAA- [^370 - kolz]^l\p exp [i{kx)] . 



(27r)4 J {kl - kl) fc2 + K2 
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After this is contracted with the unit matrix we get for the corresponding contribution into 
the right-hand side of the equation to be written in place of Eq. the expression 

_J_[ exp [ijkx)] , 

Once the Hermite functions (fT7|) restrict x_[_ in integrals like ()29|) and ((221) ^^e region 
inside the Larmour radius, the region fc^ ^ = m?B /Bq in the integral is important. 
There, however, ^2 disappears due to the exponential factor in (jHH|) and we are left with the 
contribution, the same as the one coming from the free photon propagator. Moreover, as 
the light-cone singularity is formed exclusively due to integration over near-infinite values of 
all the four photon momentum components, Eq. (PT|) behaves like the same as ()24p. 

near the light cone = 0. 

We need, however, to also estimate the contribution immediately close to this singularity. 
To this end, let us disregard the spatial dispersion of the dielectric constant in the transverse 
plane, i.e. take k,2 at the value fc^ = in (jnTjl . By doing so we essentially underestimate 
the contribution of the mode-2 photon as a carrier of the electromagnetic interaction into 
the attraction force between the electron and positron near the light cone, because we keep 
the term linearly growing with the field in the denominator for large k_i_, where it in fact 
disappears. This approximation does not affect the light-cone singularity, which remains 
l/x^, but makes the screening correction to the singular part larger than it is. We shall 
see, nevertheless, that even within this approximation, with the screening overestimated, 
the effect of the latter is small. Once our working domain is restricted to the intervals 

— t"^ much closer to the light cone than the Compton length, we may keep to the condition 
1^0 ~ ^3 1 integral ()91|) . Then k,2 (IHSl) should be taken in ()91|) as 

2aBm? 

= r—. (92) 

-kBq 

Then fl9H) becomes the well-known expression for the free propagator of a massive particle 
with the mass squared = — ^2. (To avoid a possible misunderstanding, stress that 
this mass should not be referred to as an effective photon mass. The mass of the photon 
defined as its rest energy is always equal to zero, corresponding to the fact that the point 
^0 = ^3 = ^1 = ^2 = is solution to the dispersion equations fc^ + na{k) = 0. The "mass" 
M appears in the ultraviolet, and not infrared regime.) In the adiabatic approximation 
of Subsection C of Section II the dimensional reduction yields again the prescription to 
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disregard the dependence on x± in it by setting x_l = 0, = —s^. Then the contribution 
from (|^T| is 



/Ci (Ms 



, Ms 
In 



47r2s -^^-^"^ 4vr2s2 2 ^ ' ' ^^^^ 

Here /Ci is the McDonald function of order one, and we have pointed its asymptotic behavior 
near the point = — t^ = +0 . According to (f73j) near the lower edge of the normalization 
volume and with the magnetic field ()76p the quantity sM = (2a/7r)^/^ makes 0.068, and 
hence the second term inside the brackets in (jHSj) is only —7.8 x 10~^. Therefore the screening 
effect, although overestimated, is still negligible, the contribution of D2 making one half of 
the full contribution of the free photon propagator considered above. The one half originates 
from the absence of the factor 2 that appeared above when [7i]AA<= [7«]AiAP in later 
contracted with the unity / to lead to (jUUj) : "^i^Q ^ gulili = 2. The other half comes from 
the contribution of Di and D4. 

The quantity D3 contains only k± components that give rise to [fc*7j]AA<= [^"'TjlA'Ap, hj = 
1,2, in an equation to appear in place of (jHSj), and consequently contribute only to the 
nondiagonal in Landau quantum numbers part of the Bethe-Salpeter equation like (001) , 
that does not survive in the limit of high magnetic field. On the contrary, the contributions 
of Di and D4 go to the diagonal part. This occurs because these contain the components 
fco and fcs carrying the matrices 70 and 73 that may lead to the term diagonal in Landau 
quantum numbers, as explained when passing from to ()4fi|l and ()47j) . It follows from 
(jHTjl . (jHH|) that the common contribution from Di and D4 in the (0,3) subspace is 

determined by the expression 

{bwy + F " ki-kr ^^^^ 

Then the counterpart of (j90j) reads 

_ 1 f [kplo - k^-fslxx- [fco7o - hl3]^i\p exp [ijkx]] . 

{27rrJ {kl-kl) k^ + K, ^^'^ 

and the counterpart of (j9T| becomes (again with the disregard of the spatial dispersion 

across the magnetic field already done when writing Eq. ()85|l ) 

1 [^M^^^^i 



A{2tiY J fc2 ^i47r2x2 ' 

where 
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in view of ()85p. For the fields as large as -B = B^it fl76p the number A is very close to 
unity: ^ = 1 — 0.04. (Its difference from unity is the measure of the anti-screening effect of 
the running coupling constant a/ A for large magnetic field due to the lack of asymptotic 
freedom in pure quantum electrodynamics). 

We conclude that the vacuum polarization does not any essentially affect the falling to 
the center and hence the estimate of the ultimate magnetic field. This contradicts the 
prescription to replace a a/2 in the expression for the latter that would result if one 
applied the corresponding conclusion from Ref. [24 1 to the problem under consideration. 
The point is that in Ref. the contribution of D2 is completely disregarded for the reason 
that the term (j92j) . linearly growing with the magnetic field, is in the denominator of D2. 
We saw above that that this cannot be done: it essentially contributes to the falling to 
the center asymptotic regime of s ^ lO^^^m^^, where the probability to find the system is 
concentrated. 

Gathering the results of the present consideration together we conclude that the effect 
of the vacuum polarization leads, in the approximation where the spatial dispersion in the 
orthogonal direction is neglected, to the replacement of Eq. (j58|) by the following two- 
dimensional Bethe-Salpeter equation for high magnetic field limit including an external 
electric field and the effects of the vacuum polarization 

4'-eii|(x|j)-m]^^ [i4" + eii,(xj)-mj^^ [Xo,pr,o,p^Ml,x\\)U 

^\\l^,u 1 



□9 , _PL 

^ + (eB)2 




p2 

^2 I _Sl +2 



+ 



[eB) 



[Xo,p-o,pp(a;|j,a:J)]^^,(9^ 



Here the action of the derivatives over t and z does not extend beyond the braces, id\\ = 
i7o<9f + i73<92, □2 = 5^ — df. Remind that t = Xq — Xq, z = x^ — x^. For no electric field case 
the equation that follows from (PH|) for the singlet component to substitute for (jHOl) is 

(-□2 + m^)$(M) = ^ L + ^ , ^ U(M)(99) 

Finally, the Bessel equation for the (1,1) rotationally invariant solution now becomes 

d2$ ld$ , 2« /I .A . . 

-— + m2$ = — - + M/Ci(Ms) <!>. (100) 



ds^ s ds TVS \ s 
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We neglected the difference of A from unity. 



2. Mass corrections 

Mass radiative corrections sliould be taken into account by inserting tlie mass operator 
into tlie Dirac differential operators in the l.-h. sides of the Bethe-Salpeter equation ^ or 
()47|1 . We shall estimate now, whether this may affect the above conclusions concerning the 
positronium mass compensation by the mass defect. 

In strong magnetic field the one-loop calculation of the electron mass operator leads to 
the so-called double-logarithm mass correction growing with the field B as 3^ 

m = m ( 1 + — In^ ^ ) . (101) 



V 47r Bo/ 

For B ~ i?uit the corrected mass makes m = 3.45m. This implies that the mass annihilation 
due to the falling to the center is opposed by the radiative corrections and requires a field 
somewhat larger than ()7(i|l . To determine its value, substitute m (jlOlj) for m and Lb = 
(ei?)~^/^ for So into equation (f?^ with n = 1. The resulting equation for the ultimate 
magnetic field, modified by the mass radiative corrections, -Bcorr, 

B B. 




has the numerical solution: -Bcorr — 13 -Buit- 

When going beyond the one-loop approximation by summing the rainbow diagrams two 
different expressions for m were obtained by different authors. Ref. reports 

_ f ct 2 B \ 

m = mexp — In — . (103) 
V47r Bo J 

The use of this formula analogous to the above gives rise to an increase of the ultimate value 
by two orders of magnitude: -Bcorr = 3.5 x 10^i?uit, whereas the use of the result of Ref. 

in 

m = y ^ (104) 

cos ( ^/^In -^ 



2n Bq ^ 

would leave the ultimate value practically unchanged: -Bcorr = 1-5-Buit- Finally, if the vacuum 
polarization is taken into account while summing the leading contributions to the large field 
asymptotic behavior of the mass operator, the following result 

in 

m = 7 r 5- (105) 

1-f (in^-^E)lnf 
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was obtained, from where the double logarithm is absent due to the effect of the term ()92|) 
in the photon propagator when substituted into electron-photon loops. The use of ()105|) 
would result in: -Bcorr = 3-Buit- 

Anyway, we see that the mass correction, increasing the ultimate value i?uit by at the 
most two orders of magnitude, is not essential bearing in mind the huge values (|76|) of the 
latter. Moreover, basing on the most recent results concerning the mass correction [s^ we 
conclude that the latter do not affect the value of the hypercritival field obtained above ()76|) 
practically at all. 



IV. SUMMARY AND DISCUSSION 



In the paper we have considered the system of two charged relativistic particles - especially 
the electron and positron - in interaction with each other, when placed in a strong constant 
and homogeneous magnetic field B. The Bethe-Salpeter equation in the ladder approxima- 
tion in the Feynman gauge is used without exploiting any non-relativistic assumption . We 
have derived the ultimate two-dimensional form of the Bethe-Salpeter equation, when the 
magnetic field tends to infinity, with the help of expansion over the complete set of Ritus 



matrix eigenfunctions 19|. The latter accumulate the spacial and spinor dependence on 



the transversal-to-the-field degree of freedom. The Fourier-Ritus transform of the Bethe- 
Salpeter amplitude obeys an infinite chain of coupled differential equations that decouple 
in the limit of large B, so that we are left with one closed equation for the amplitude com- 
ponent with the Landau quantum numbers of the electron and positron both equal to zero, 
while the components with other values of Landau quantum numbers vanish in this limit. 
The resulting equation is a differential equation with respect to two variables that are the 
differences of the particle coordinates: along the time t = Xq — Xq and along the magnetic 
field It contains only two Dirac matrices 70 and 73 and can be alternatively 

written using 2x2 Pauli matrices. The term responsible for interaction with a moderate 
electric field E directed along B, ^ _B, is also included and does not lay obstacles to 
the dimensional reduction. By introducing different masses the resulting two-dimensional 
equation may be easily modified to cover also the case of an one-electron atom in strong 
magnetic field and/or other pairs of charged particles. 

It is worth noting that the two- dimensionality holds only with respect to the degrees of 
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freedom of charged particles, while the photons remain 4-dimensional in the sense that the 
singularity of the photon propagator is determined by the inverse d'Alamber operator in the 
4-dimensional, and not two-dimensional Minkowsky space. (Otherwise it would be weaker). 

We have made sure that in the case under consideration the critical value of the coupling 
constant is zero, = 0, i.e., the falling to the center caused by the ultraviolet singularity 
of the photon propagator as a carrier of the interaction is present already for its genuine 
value a = 1/137, in contrast to the no-magnetic-field case, where > 1/137. If the 
magnetic field is large, but finite, the dimensional reduction holds everywhere except a 
small neighborhood of the singular point s = 0, wherein the mutual interaction between 
the particles dominates over their interaction with the magnetic field. The dimensionality 
of the space-time in this neighborhood remains to be 4, and its size is determined by the 
Larmour radius Lb = {eB)~^^'^ that is zero in the limit B = oo. The latter supplies the 
singular problem with a regularizing length. The larger the magnetic field, the smaller the 
regularizing length, and the deeper the level. 

We have found the ultimate magnetic field that provides the full compensation of the 
positronium rest mass by the binding energy, and the wave function of the corresponding 
state as a solution to the Bethe-Salpeter equation. This state is described in terms of the 
theory of the falhng to the center, developed in Q , as a " confined" state, different from 
the usual bound state. The appeal to this theory is necessitated by the fact that the falling 
to the center draws the electron and positron so close together that the mutual field is so 
large that the standard treatment may become inadequate. The ultimate value is estimated 
to be unaffected by the radiative corrections modifying the mass and polarization operators. 

In spite of the huge value, expected to be present, perhaps, only in superconducting cosmic 
strings [3| , the magnetic field magnitude obtained may be important as setting the limits of 
applicability of QED or presenting the ultimate value of the magnetic field admissible within 
pure QED. The point is that at this field the energy gap separating the electron-positron 
system from the vacuum disappears. An exceeding of the ultimate magnetic field would 
cause restructuring of the vacuum. The question about the vacuum restructuring typical of 
other problems - with or without the magnetic field, where the falling-to-the-center takes 
place: the supercharged nucleus js^, Q and a moderately charged nucleus with strong 
magnetic field js^, is discussed in the two adjacent papers jly, Q|- The formal mechanisms 
that realize the magnetic field instability and may lead to prevention of its further growth via 
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the decay of the " confined" state found here require a further study and will be considered 
elsewhere. 
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